Force Dynamics in Weakly Vibrated Granular Packings 
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The oscillatory force F£ c on the bottom of a rigid, vertically vibrated, grain filled column, reveals 
rich granular dynamics, even when the peak acceleration of the vibrations is significantly less than 
the gravitational acceleration at the earth's surface. For loose packings or high frequencies, F^ c, s 
dynamics are dominated by grain motion. For moderate driving conditions in more compact samples, 
grain motion is virtually absent, but F^ c nevertheless exhibits strongly nonlinear and hysteretic 
behavior, evidencing a granular regime dominated by nontrivial force-network dynamics. 

PACS numbers: 45.70.-n,05.45.-a 



Granular media consist of macroscopic solid grains 
which interact via dissipative, repulsive contact forces. 
Thermal energy is inconsequential, and granulates jam- 
in random configurations unless sufficient mechanical en- 
ergy is supplied, for example by shearing or shaking 
For sinusoidally, vertically vibrated granular media the 
driving strength is characterized by the nondimcnsional 
acceleration amplitude, T = A(2irf) 2 /g, where A is the 
displacement amplitude, / is the oscillation frequency, 
and g is the gravitational acceleration at the earth's sur- 
face. In the well-studied case of T > 1, grains periodically 
lose contact with and subsequently impact the oscillat- 
ing container. The collision, as well as the accompanying 
relative motion between grains and wall, injects energy 
into the system that drives grain rearrangements 0, U 

In this Letter we explore the T < 1 regime by measur- 
ing the oscillatory force on the container bottom, F£ c , 
in a weakly vibrated column filled with grains of total 
mass M [Fig. 1(a)]. F£ c exhibits rich dynamics. Un- 
der two conditions the dynamics is dominated by grain 
rearrangements: (i) Loose packings may compact when 
they are first subjected to weak vibrations, which leads 
to intermittent burst in F£ c . (ii) For large driving fre- 
quencies, grains slide periodically with respect to the con- 
tainer even at low T (e.g. V « 0.1 at / = 900 Hz), which 
leads to strongly nonlinear behavior in F£ c . 

With moderate driving in compacted samples, how- 
ever, relative grain motion is minute but variations in 
the force configuration remain substantial as indicated 
by nontrivial changes in F b ac . Weak vibrations thus ex- 
cite strongly nonlinear and glassy dynamics of the force 
network. Such force variations are allowed, since for a 
terrestrial granulate the grain and deformation scales are 
separated by many orders of magnitude Q : a 700 /jm di- 
ameter bronze sphere is compressed only ^100 nm under 
the weight of 1000 additional identical spheres 

Figures ^b,c) illustrate the main features of the 
first harmonic of F£ c divided by T [see Eq. QJ], f° r a full 
column of grains undergoing slow triangular sweeps of T. 
For a solid mass placed on the bottom plate, F\ is inde- 
pendent of r. For grains, however, F\ depends strongly 
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FIG. 1: (a) Schematic of the experiment showing the piezo 
force sensor mounted between the "bottom" and the "base" 
(diameter x height: 32 x 12 mm and 89 x 62 mm respectively) , 
with a cylindrical tube (inner/outer diameter 30/55 mm, 
height 113 mm) which is attached only to the base. The small 
~ 100 /im gap between the tube and bottom plate prevents 
grains from becoming trapped. Shaded parts are rigidly con- 
nected and move in unison. The sensor signal is used to obtain 
F§ c , the vertical AC force exerted on the bottom plate by the 
grains. (b,c) Nonlinear response of a column filled with 200 g 
of 0.61 — 0.70 mm diameter bronze particles under sweeps 
of the vibration amplitude F for / = 80 Hz; F\ denotes the 
calibrated ratio of the first harmonic of F£ c to T (Eq. ^ . 

and nonlincarly on T [e.g. Fx(T = 0.5) « 2Fi(r = 0.05)], 
even for low driving frequencies. The spikes in F\ during 
the initial ramp are caused by compaction of the material. 
The asymmetry of -Fi(r) indicates hysteresis and mem- 
ory in compacted samples. We stress that the strength 
of these features does not vary significantly for driving 
frequencies from 16 to 300 Hz. The phenomena evident 
in Fig. ^b) are essentially quasistatic and are not asso- 
ciated with the excitation of sound waves (see below). 

Experimental Setup — Nearly-spherical bronze parti- 
cles sieved between 0.61 and 0.70 mm are poured into a 
smooth cylindrical tube with a detached bottom which 
is supported by a rigid piezo-electric force sensor (stiff- 
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ness 2.5 GN/m). The entire assembly is vertically oscil- 
lated with a small sinusoidal displacement [Fig. 1(a)]. An 
accelerometer attached to the tube measures the time- 
resolved acceleration 7(f) which, for most driving condi- 
tions, is harmonic, equaling rsin(27r/t). The measured 
force is sensitive to temperature drift. The entire as- 
sembly is therefore placed in a temperature controlled 
enclosure maintained slightly above room temperature 
(temperature fluctuations ± 10 mK, humidity 5-10%); 
grains are equilibrated in the enclosure prior to use 6]. 

The deflection of the relatively compliant force sensors 
used in most previous studies is large compared to the 
deformation of hard grains like steel or glass; the gran- 
ular force configuration is then completely altered due 
to relative motion between force probe and grains. In 
contrast, our tube/sensor assembly is effectively a solid 
container since the maximal deflection of the piezo is less 
than 1 nm, which ensures that the measured force varia- 
tions are intrinsic to the granular medium 0, S @ ■ 

In our experiment, the force sensor measures the total 
AC force, F^ tal , which is the sum of the inertial force 
generated by the acceleration of the bottom plate and 
sensor, with effective mass Mo, and the AC bottom force 
F b ac resulting from the acceleration of material in the 
column (which can have harmonics): 

00 

F t a c tal = TM g sin(2^/t)+r £ F n sm(2nnft-cb n ) . (1) 

n=l 

To calibrate the signal, F^ tal was measured for a range 
of / and r, both with and without solid test masses at- 
tached to the bottom plate. The value of F^ tal for the 
empty system allows us to subtract the term oc TMog, af- 
ter which F\ is found to be proportional to the test mass 
and independent of T (the definition of F n isolates the 
trivial scaling with T). The higher harmonics F2, F3, . . . 
are negligible in this case. 

We now consider the AC vertical components of the 
frictional wall-force F w and the bottom force F^ in a grain 
filled column. We will distinguish a contact regime where 
grains do not slide with respect to the column, and an 
impact regime where they do. In the impact regime, 

MgTsm(2nft) = F^ + F b ac , (2) 

and our measurements probe how the vibration induced 
variations of the grain weight, M gT sm(2% ft), are dis- 
tributed between wall and bottom. However, Eq. J5J is 
violated in the impact regime. We have checked the valid- 
ity of Eq. J2J by placing a separate, closed bottom, grain 
filled container directly on the bottom plate ^(J > so that 
the sensor measures the sum F^ c + F b ac ; when the sensor 
signal remains purely harmonic, the system is in the con- 
tact regime, when (strong) nonlinearities are detected, it 
is in the impact regime. As we will show for the contact 
regime, F b ac grows nonlinear ly with grain mass and T, 
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FIG. 2: Spikes observed for / = 80 Hz and M = 200 g in an 
initially low density packing. Sweeps of increasing magnitude 
in F illustrate that spikes only occur in "fresh" territory. 

exhibits hysteresis and memory effects, and higher har- 
monics play a role — very different behavior than for a 
solid mass. However, we first discuss grain motion dom- 
inated compaction and impact. 

Compaction - - "Spikes" , such as those shown in 
Fig- 12b), occur when T is ramped up in loosely packed 
samples, which are formed by placing the end of a funnel 
on the bottom of the container, filling the funnel with 
material, and then slowly retracting the funnel. Figure[21 
illustrates that spikes only occur in "fresh" territory, i.e., 
when r is increased beyond its previous maximum value. 
After r has been swept up to a value near one, we refer 
to the sample as fully annealed. 

During a spike, which typically lasts for 1000's of os- 
cillation cycles, <pi shifts significantly (indicating dissi- 
pation) and F b ac is strongly non-sinusoidal. The specific 
r values where spikes occur vary from run to run, and 
hence are not resonant effects. Spikes are apparently due 
to compaction, since the free surface is lower after a spike 
occurs, and gently poured columns with lower initial den- 
sity produce more spikes than less gently prepared ones 
with higher initial density. Spikes and their dependence 
on r apparently result from the full mobilization of fric- 
tional forces at the walls in granular columns 0, EI ■ Un- 
der these conditions, vertical frictional forces are near 
their maximal values so that vibrations can cause (mi- 
cro) slippage, and stronger vibrations may cause further 
compaction. The behavior of spikes is not frequency de- 
pendent: qualitatively similar spikes occur for / between 
16 and 300 Hz, and switching to a different frequency in 
a partially annealed sample does not appreciably alter 
the r value separating fresh and annealed states [6j. 

Impact vs Contact — For compacted samples we dis- 
tinguish between the contact regime where Eq.[3is satis- 
fied, and the impact regime where Eq.[5]is violated. Fig- 
ure a) illustrates that for low frequencies F\ increases 
smoothly with T, while for higher frequencies there is a 
sudden upturn and a peak. In the vicinity of and above 
this transition, j(t) and F b ac are strongly anharmonic — 
as when grains periodically collide with the bottom for 
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FIG. 3: Contact and impact regimes (M — 200 g): (a) 
Fi(T,f) is smooth at low frequencies, but increases abruptly 
and peaks for higher frequencies (data is vertically offset for 
clarity) . (b) Greyscale intensity plot of nonlinearity Af(T, f) = 
Fi(r,/) - Fi(0,/) (white, A/" = 0; black, JV positive) indicat- 
ing the contact (C) and impact (J) regimes. Dots indicate 
the initial sharp increase in F\. 



FIG. 4: Nonlinearity of F\ with F and total grain mass M. 
(a) Fi(T)/Mg for various filling fractions (/ = 80 Hz). For 
larger M and fixed F, the weight fraction on the bottom 
decreases, (b) Fi as function of M in the linear regime 
(r = 0.05) for / = 16, 80 and 300 Hz, compared to a 
linear response (straight line) and a Janssen-like response 
Fi/F sat = 1 - exp{-Mg/F sat ) for F sat = 0.435 N. 



r > 1; this is the impact regime. Below the transition, 
j(t) remains sinusoidal; this is the contact regime. 

To check that Eq. J2J is valid in the contact regime 
(i.e. no slipping) but not in the impact regime, we re- 
placed the open-bottom tube with a similarly sized, but 
closed-bottom grain filled tube mounted directly on the 
bottom plate. In this configuration, F b ac is insensitive 
to smooth force transfers between the wall and the con- 
tainer bottom which allows Eq. to be checked directly. 
In the impact regime F£° and "f(t) are strongly anhar- 
monic in violation of Eq. ; in the contact regime they 
remain sinusoidal with F\ essentially independent of F in 
agreement with Eq. |J2J. The physics underlying the force 
dynamics in the contact regime is thus a smooth, periodic 
transfer of grain weight between wall and bottom. 

Figure displays the strength of the nonlinearity 
of F\ as function of Y and / — the rapid increase in non- 
linearity marks the onset of impact. Note that at high 
frequency kHz) , the impact regime occurs for sur- 

prisingly small r (wO.l). Apparently, impact for T < 1 is 
due to the excitation of resonant granular sound waves. 
Typical sound speeds are of the order of 100 ms _1 so 
in our 10 cm deep column we expect a resonant response 
around 1 kHz. This picture is consistent with findings of 
Yanagida et al. in studies of the resonant response of 
grain filled closed bottom containers for small T, and it 
is also consistent with the shift of the impact transition 
to higher frequencies for smaller M (this also excludes a 
trivial resonance of the apparatus 6j). The T dependence 
of this transition is not fully understood. 

Contact Regime Nonlinearity — We now explore the 
nonlinear response of fully annealed samples in the con- 
tact regime as a function of M and T. The frequency 
is fixed at 80 Hz, since in the contact regime the grain 
response varies only weakly with /. 

Figure |3J a) illustrates that the small mass behavior is 
independent of T, and that the grains are supported en- 
tirely by the bottom since F\ w Mg. For larger masses, 



wall forces start to play a role since F\ < Mg and the 
material's response becomes increasingly nonlinear with 
T: the exchange between bottom and wall forces under- 
lies the nonlinearity of i*i(r) for large M, This strong 
nonlinearity is accompanied by a modest increase in the 
harmonic distortion 01 of F b ac and a similarly modest 
shift in For example for T = 0.7 and M — 200 g, 
F\ is doubled from its low Y value, while the harmonic 
distortion is only 6% and <f>\ ~ 6°. Underlying the har- 
monic distortion is the nonlinear (Hertzian) contact-force 
law which makes grain contacts stiffcr when under larger 
pressure 0; the waveform distortions we have obtained 
for F£ c (t) are consistent with this picture @]. 

For small T the response is linear in T: for T as 0.05, 
F fc ac is harmonic (< 1% distortion) and in-phase with the 
acceleration, and Fi varies less than 1% for < T < 0.1. 
We study Fi(M) as a function of increasing mass M at 
r« 0.05 by incrementally pouring grains from a height of 
approximately 10 cm above the grain surface. Figure|l| b) 
illustrates that F± grows proportionally with M for small 
masses, but then rapidly saturates to F( at ~ 0.435 N. 
F\(M) is only weakly frequency dependent, again indi- 
cating that in the contact regime a well-defined quasi- 
static regime is probed. Note that for all / a small over- 
shoot occurs for intermediate values of M w 100 g. 

This behavior is reminiscent of the Janssen effect for 
which the DC bottom force F b dc goes as HEUQ: 

Ft = [1 - exp (-Mg/Ft t )} , (3) 

where Fjj£ t is the saturation force. Figure Elb) shows, 
however, that the AC bottom force F\ significantly de- 
viates from the static Janssen result. We conclude that, 
even in the limit of weak vibrations, the force variation 
F£ c is not simply related to the steady force F b dc . 

Hysteresis and Memory — So far, our data is consis- 
tent with a simple nonlinear response of the material, 
but Fig. |SJ a) indicates that when T is ramped up and 
down in the nonlinear regime, F£ c is also hysteretic. The 
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FIG. 5: Memory effects (M = 200 g and / = 80 Hz), (a) 
Hysteresis loop, (b) Hysteresis persists when T is suddenly 
set to zero and then rapidly ramped back up, indicated by 
the circle (see inset for details of ramp), (c-f) Two examples 
of subtle memory effects. 

magnitude of hysteresis is only weakly dependent on the 
driving frequency and increases with M and T similar to 
the magnitude of overall nonlinearity. The hysteresis is 
nearly independent of the sweep speed (for sweep dura- 
tions longer than ss 100//), and the force configuration 
thus depends on the driving history. Fig.[SJb) illustrates 
that these configurations can be "frozen" , since the sys- 
tem returns to the upper branch of the hysteresis loop 
after the driving is switched off and then rapidly ramped 
up again. The system exhibits additional subtle memory 
effects: after a fully annealed system is subject to a num- 
ber of small amplitude sweeps in T, and then T is ramped 
beyond the peak value of the small sweeps, a clear "kink" 
in the iq curve is exhibited [Fig. E2c,d)]. When a fully 
annealed system is driven at a fixed T, it "remembers" 
this value when T is rapidly decreased and then ramped 
past the initial fixed value [Fig.Efe,f)]. 

Discussion — Our experiments exhibit rich dynam- 
ical behavior of weakly excited (r < 1) granular media 
which can be dominated by either grain motion or by con- 
tact force variations. That grain motion and compaction 
occur in loose samples vibrated at low frequencies and 
r < 1 is maybe not surprising, although we are unaware 
of systematic studies of compaction in this regime [3j . In 
compacted samples, both grain dominated (impact) and 
force dominated (contact) regimes can be distinguished 
with / playing a crucial role. For low /, the system is in 
the contact regime, and we observe little systematic vari- 
ations with /. We therefore interpret the nonlinearities, 
memory effects, and hysteresis of the forces in response 
to variations in the driving strength as quasistatic, and 
not related to the excitation of sound waves. We have 
not identified any theoretical or numerical descriptions 



of these surprisingly strong effects. Exploratory experi- 
ments in various columns with rough walls and for parti- 
cles of different sizes further illustrate the robustness of 
these phenomena Q. 

How should weakly vibrated granular systems be 
viewed? The contact laws between elastic bodies display 
a variety of nonlinear and hysteretic behaviors |a, ll 51 ] . A 
weakly driven granular assembly apparently amplifies the 
local nonlinearities present in the hertzian contacts and 
the friction law. We propose that force networks "acti- 
vated" by weak vibrations explore many different config- 
urations consistent with the overall boundary conditions 
for the stress 0, 0] . Such activated force networks could 
possibly play a role in creep flows, which occur far away 
from shear zones, and more generally in any granular sys- 
tem in which tiny relative grain motions are excited. In 
this sense, weakly driven granulates cannot be thought 
of as ordinary solids. 
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